SCIENTIFIC 

REPORTS 



ul7 




OPEN 



SUBJECT AREAS: 
COEVOLUTION 
BEHAVIOURAL ECOLOGY 
SOCIAL EVOLUTION 

Received 
21 January 2014 

Accepted 
13 February 2014 

Published 
4 March 2014 



Correspondence and 
requests for materials 
should be addressed to 
M.A.C. (abouchakra@ 
evolbio.mpg.de) 



Plastic behaviors in hosts promote the 
emergence of retaliatory parasites 



Maria Abou Chakra\ Christian Hilbe^'^ & Arne Traulsen^ 



^ Evolutionary Theory Group, Max Planck Institute, D-24306 Plon, Germany, ^Program for Evolutionary Dynamics, Harvard 
University, Cambridge, MA 02 1 38, USA. 

Mafia like behavior, where individuals cooperate under the threat of punishment, occurs not only in 
humans, but is also observed in several animal species. Observations suggest that avian hosts tend to accept a 
certain degree of parasitism in order to avoid retaliating punishment from the brood parasite. To 
understand under which conditions it will be beneficial for a host to cooperate, we model the interaction 
between hosts and parasites as an evolutionary game. In our model, the host's behavior is plastic, and thus, 
its response depends on the previous interactions with the parasite. We find that such learned behavior in 
turn is crucial for the evolution of retaliating parasites. The abundance of this kind of mafia behavior 
oscillates in time and does not settle to an equilibrium. Our results suggest that retaliation is a mechanism 
for the parasite to evade specialization and to induce acceptance by the host. 

Host-parasite coevolution involves constant adaptation and counter adaptation; the host acquires defenses 
to deflect the costs imposed by the parasite; as a response, the parasite evolves new ways to escape from 
host defense. Parasites are considered as 'free-riders' since they exploit the host in order to increase their 
own reproductive success. For example, brood parasites will lay their eggs in the host's nest and evade all parental 
care. This host-parasite interaction causes a cascade of co- evolutionary traits such as mimicry, signatures, and 
timing defenses^"^. However, in spite of their abilities to defend themselves, some hosts are observed to tolerate the 
parasite. Brood parasites and their hosts are a model system, but despite ample experimental evidence describing 
this CO -evolutionary process, it is still not clear under which conditions the hosts are compelled to accept 
parasitism'*"^^. 

Cooperating with a parasite has detrimental consequences. Two hypotheses have been proposed to explain 
when accepting behavior is evolutionarily possible^'^^. The evolutionary lag hypothesis posits that hosts would 
naturally evolve to reject the parasites. However, there is a lag in their response and insufficient amount of time 
has passed for the adaptation to occur. The equilibrium hypothesis posits that the costs associated with rejection 
outweigh acceptance. Rejection is costly, and thus, less adaptive. Both hypotheses have some support, but 
conclusive evidence is wanting^"^^. 

Alternatively, Zahavi proposed a 'mafia hypothesis' which posits that hosts accept because they fear retaliation 
by the parasite^\ Empirical evidence shows that both great spotted cuckoos {Clamator glandarius) and brown- 
headed cowbirds {Molothrus ater) exhibit such retaliatory behavior. These parasites have different strategies and 
host types, yet, both were observed to punish their respective hosts {e.g. magpies. Pica pica, and prothonotary 
warblers, Protonotaria citrea) for rejecting their eggs by destroying entire successive nests. In order to avoid 
repercussions from the retaliating parasitic behavior, hosts tend to accept a certain degree of parasitism^^'^^. The 
mafia hypothesis works when retaliation is sufficiently costly to the host in comparison to the consequences of 
accepting parasitism. In particular, this assumes that the host's young are raised alongside of the parasite's young, 
which is the case in these examples However, the behavioral consequences of retaliatory punishment are still 
unclear. Here, we investigate the emergence of retaliatory behavior and explore the host's response to such 'mafia- 
like' behavior in parasites by an analytical model and numerical simulations. 

Direct behavioral adaptations and counter-adaptations render the brood host-parasite system described above 
ideal for theoretical investigations^'^^"^^. Robert et al. showed, by including repeated interactions, that parasite 
retaliation oscillated temporally; they observed that retaliation evolves when host rejection is abundant^^. 
However, typical behavioral models^^"^^ have only considered two types of hosts (rejecters or accepters) and 
parasites (retaliate rs or non-retaliators). In these previous models, host do not exhibit any plastic response to 
parasites. Since hosts are incapable of adapting their behaviors between clutches or after an interaction with a 
parasite, previous models could not address whether the retaliatory mafia truly influences the decision of the host 
in subsequent interactions. In contrast, in nature some hosts change their behaviors by learning within a single 
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breeding season^'^^'^^. Hence, in our model, we allow hosts to react on 
their past experience with parasites: A host could accept parasitism 
unconditionally (first clutch) or conditionally (reject in first clutch 
and accept in second clutch). Spiteful behaviors can emerge, parasites 
who retaliate, but do not return to the host. Also, opportunistic 
behaviors can emerge, parasites who return to the host without 
retaliating. Given all those behaviors interacting together, will 
mafia-like' behavior evolve? When should a host accept condition- 
ally versus unconditionally? Additionally, how much parasitism 
should a host accept? To address these questions, we construct a 
model that captures these more complex host-parasite interactions. 

Results 

The model. We model the interaction between hosts (e.g. magpies. 
Pica pica, and prothonotary warblers, Protonotaria citrea) and 
parasites (e.g. cuckoos, Clamator glandarius, and cowbirds, 
Molothrus ater) as an evolutionary game^^'^^'^^ These species tend 
to interact with each other over several seasons. In each season a 
parasite can lay more than a single egg in the host's nest without fear 
that the parasite nestling will evict its siblings out of the nest; both the 
young of the host and parasite are raised together^'^^'^^. Thus, in this 
system there are repeated interactions within a breeding season. We 
classify the host-parasite interactions into four stages: (i) During the 
first parasitism stage the parasite lays a certain number of eggs in the 
nest and the host will choose to accept these eggs or to reject them, (ii) 
If the host rejects, the parasite could depredate and destroy the host's 
nest during the first depredation stage. This stage occurs after the 
host rejects, which could be any time between egg laying and nestling 
hatching^^. (iii) The second parasitism occurs only when the parasite 
returns to the same host to lay additional eggs after rejection or after 
the host reclutches. (iv) The second depredation stage occurs after a 
host rejects the returning parasite. 

In the model, a parasite's behavior is determined by the number of 
eggs laid per host nest, whether it is a retaliator (mafia) or not, and 
whether it returns to the same host or not. Similarly, the host's 
behavior is determined by the number of parasitic eggs tolerated, 
whether retaliating parasites are accepted or not upon their return, 
whether non- retaliating returning parasites are accepted or not, and 
whether the host ejects the parasitic eggs or abandons the entire nest. 
A host's behavior can differ between stages. For instance, a host can 
reject in the first parasitism stage, but choose to accept in the second 
parasitism stage. These behaviors can become complex quickly, and 
thus, it is important to distinguish between each response that is 



possible by a parasite {e.g. number of eggs laid per nest, retaliation, 
returning to the same host) and host (e.g. acceptance, rejection, eggs 
tolerated, accept retaliators or not). We assume that a host will lay 
one clutch during the breeding season unless they are depredated on, 
with a maximum of three clutching attempts in a breeding season^^'^^. 

This is captured by an asymmetric evolutionary game played 
between one parasite and multiple hosts selected at random from 
their respective populations; this setup mimics natural scenarios 
where parasites, depending on the species, can lay between 20 and 
40 eggs per breeding season^ Using evolutionary game theory, we 
analyze host-parasite interaction at each stage. Over time, behaviors 
evolve based on the fitness^^. In each season, a parasite has the capa- 
city to lay a certain number of eggs that is allocated among hosts. 
Non -mafia parasites lay a total of f^^ eggs while mafia parasites lay a 
total of l^M eggs. When [^m > I^n the emergence of mafia behavior is 
trivial, and will not be discussed herein. Instead, we assume Pm — Pn^ 
that is mafia-behavior may be costly for the parasite (for example, 
because the parasite may mistakenly eliminate its own offspring, or 
because of cognitive and energy costs^^). A host lays bh eggs in each 
clutch, and under parasitism incurs a cost Cp if parasitism occurs, c„ if 
a nestling is raised, and if forced to re- clutch. To calculate the 
average fitness of a parasite we consider the average fraction of 
accepted eggs reared by the host, whereas to calculate the average 
fitness of a host we consider the average number of its eggs success- 
fully reared. 

For a concrete analysis we consider two models, a minimalistic 
version with few strategies that can be analyzed analytically and a 
more complex version which we analyze with simulations. Both 
models allow hosts to change their behavior between each interaction 
(plastic response). 

Minimalistic model. We first consider the case where the parasite 
lays a single egg per nest, bp= I and only consider three types of host 
and two types of parasites. Host types are accepters (hosts that accept 
after first parasitism stage), conditional accepters (hosts that accept 
only after retaliation), or rejecters (hosts that always reject). Parasite 
types are non-mafia (non- retaliators that switch to a new host after 
rejection) or mafia (retaliators that depredate the nest and re- 
parasitize the same host). The fitness of both, host and parasite, 
depend on their evolved strategies. 

We can infer the average fitness of a host from the game tree in 
Fig. 1. We suppose that in each breeding season, the host lays bh eggs 
in a clutch which can become parasitized with a single egg, bp= 1. The 
host may accept the parasitic egg incurring a parasitism cost Cp (pam- 
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Figure 1 | Game tree for the host-parasite interaction in the minimalistic model. The host lays bh eggs in a clutch which can become parasitized with a 
single egg. The hosts fitness may be reduced by Cp (when the host is parasitized), by c„ (when the parasite's nestling is raised), and by (when the host is 
forced to re-clutch). The parasite gains the accepted egg, bp. 
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sites often remove or destroy host's eggs, typically as many eggs as 
they lay in a nest^) and a nestling cost, c„ (for raising the parasite's 
young); the parasite gains the accepted egg, bp. Alternatively, the host 
may reject the parasitic egg, but risk retaliation whereby the parasite 
destroys the whole clutch. In case of non- retaliation, the host just 
incurs a parasitism cost, without any gain for the parasite. In case 
of retaliation, the host lays a second clutch at a cost, c^, which again we 
assume is parasitized. In this situation, the host may accept the second 
parasitic egg and incur a parasitism cost, nestling cost and a reclutch- 
ing cost; the parasite gains the accepted egg, bp. Alternatively, the host 
may reject the second parasitic egg, and due to the retaliatory behavior 
of the parasite, both end up with nothing (we assume retaliators 
always return and all future clutches are destroyed). 

For the analytical model, the three types of hosts, accepters A, 
conditional accepters C, and rejecters R, have respective frequencies 
denoted by y^, yc and)/^^ such that + 7c + JKi? ^ 1- Similarly, the 
two types of parasites, non-mafia N, and mafia M, have respective 
frequencies denoted by and Xm such that x^ + Xm = 1- 

To calculate the average fitness of a parasite we compute the aver- 
age acceptance rate for each parasite type. For a non-mafia strategist, 
since they lay a single egg per nest, the average number of eggs 
accepted per parasitized host is directly proportional to accepter host 
frequencies, j/^- For a mafia strategist, on the other hand, the average 
number of eggs laid per parasitized host is ly^ + 2yc + 2yji = 2— y^^, 
whereas the expected number of accepted eggs per parasitized host is 
yA + yc- Thus, the average fitness (ratio of accepted to laid eggs) of a 
parasite is given by 

'^N=yA'PN^ 

yA+yc r> 



TlM = 



2-yA 



We infer from equations (1) that when all hosts are accepters {y^ = 1) 
the fitness of non-mafia and mafia-types are jS^and f^M^ respectively, 
such that non-retaliatory parasites have the higher fitness. When, on 
the other hand, all hosts are conditional accepters (yc = 1) the non- 
mafia parasites are always rejected and every second lay by a mafia 
parasite is accepted. In such a host population, retaliation pays. 
Finally, when all hosts are rejecters {yR= 1) both types of parasites 
fail to reproduce. 

To calculate the average fitness of a host, we first calculate the 
probability with which each parasite type visits the host. A non-mafia 
parasite needs to visit f^^ different hosts (one host for each egg), 
whereas a mafia parasite requires on average PmK2 — yA) different 



hosts (which is between Pm if all hosts are accepters, and Pm^2 if not a 
single host is an accepter and the parasite returns to all of them). 
Thus, the average fitness of a host is 



TlA = bh—Cp—Cn, 

nc = bh — Cp- 



PMXM/{2-yA) 



Ur-- 



Pn^n + Pm^m/ (2 -yA) 



(2) 



Pn^n + Pm^m/ (2 -yA) 



(bh-Cp). 



In the absence of retaliatory parasites (i.e., Xm = 0), rejecters and 
conditional accepters have a higher fitness than accepters. This rela- 
tion, however, is reversed when retaliation is common (i.e., Xm= !)• 
Moreover, conditional accepters dominate rejecters (in the sense that 
Uc ^ tir, and ttc > tir if the parasite population contains any mafia 
strategists, Xm> 0). Thus, for modeling the evolutionary dynamics, 
we simplify our model by neglecting unconditional rejecters R in the 
host population in the following. 

We model the dynamics in this host-parasite population using the 
replicator dynamics for asymmetric games^^; that is, we assume that 
the evolution of strategies can be described with the following two 
equations: 



Xm=Xm{'^—Xm){t^m — t^n), 

yA=yA{^-yA){nA-nc). 



(3) 



The dynamics of this host-parasite interaction is cyclic (Fig. 2). When 
the non-mafia parasites are in the majority, conditional accepter 
hosts have the higher fitness. As the frequency of these conditional 
accepters increases, it becomes beneficial for the parasites to retaliate 
against these C type hosts, and thus, mafia parasites increase in 
frequency. As soon as mafia parasites are common, it is optimal 
for hosts to give in without delay, leading to an increase in accepter 
hosts. This, in turn, makes it needless for parasites to retaliate, leading 
the parasite population back to the non-mafia strategy. 

All of the pure population states are unstable fixed points. 
Moreover, the cyclic dynamics has a fixed point (x^, y*^) in the 
interior of the state space. 
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Figure 2 | Evolutionary cycles in the replicator dynamics for the host-parasite interaction. The left graph shows the state space, with the horizontal axis 
corresponding to the parasite frequencies (mafia behavior is frequent at the left side, and rare at the right side), whereas the vertical axis gives the host 
frequencies (accepters are frequent at the top, and conditional accepters at the bottom). Dashed lines show the isoclines of the system, their intersection 
gives the unique, but unstable equilibrium in the interior. Even if the initial population is close to the equilibrium, evolutionary orbits follow 
periodic cycles with increasing amplitude, as shown in the right graph (parameters ^/j = 5, = Cp = c„ = 1, Pm =18 and Pn^ 20). 
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Figure 3 | Game tree of host-parasite interaction in the complex model. A host lays bh eggs in a clutch which can become parasitized by bp parasitic eggs. 
Depending on the host's threshold trait t, the host may accept the parasitism (if bp < t) incurring a parasitism cost Cp, a nestling cost c„. When hosts reject 
parasitism (if bp>t)y they risk being forced to re-clutch, resulting in costs c^. The parasite gains the accepted egg, bp. 



For the parasites, the equilibrium fraction of mafia strategists 
depends on their relative baseline fitness Pm^Pn^ and the costs 
incurred by the host, the nestling cost c„ and the re-clutching cost 
Cr. Higher re-clutching costs decrease the equilibrium abundance of 
mafia parasites. As Cr is increased, it becomes riskier for hosts to reject 
the parasite, and thus, the mere threat of retaliation is often sufficient. 
Surprisingly, the equilibrium frequency of mafia parasites increases if 
their relative baseline fitness Pm^Pn is reduced. This counterintuitive 
effect of the relative baseline fitness can be understood by analyzing 
the equilibrium frequency of accepting hosts. A low Pm reduces the 
probability that a given parasitic egg is laid by a mafia-parent, thereby 
making hosts more inclined to reject foreign eggs in their first clutch. 
This in turn promotes the spread of mafia behaviors that enforce the 
host's acceptance in subsequent clutches. Overall, this effect out- 
weighs the initial disadvantage of a lower baseline fitness Pm- 

Notably, the host's equilibrium abundance is independent of the 
costs Cp, Cr and c„. An increase of these costs therefore does not affect 
the eventual equilibrium fraction of accepters and conditional accep- 
ters in the population. Instead, changes in theses costs are only 
reflected by strategy changes of the parasites (this somewhat counter- 
intuitive effect is characteristic for mixed Nash equilibria, see for 
example^^). 

However, the interior equilibrium is evolutionarily unstable and a 
numerical analysis suggests that there is a limit cycle (as this cycle is 
close to the boundary of the system in our case, in a finite population 
extinction of one type in each population is a likely outcome). These 



dynamical results deviate from the equilibrium hypothesis: although 
a cost-benefit equilibrium exists, any small perturbation will drive 
the population away, and once left, there is no evolutionary force 
leading back. Interestingly, our model suggests that the resulting 
evolutionary lag does not diminish over time. 

Complex model. The aforementioned results give a proof of 
principle; they show how mafia behavior can emerge in a 
comparably simple setup. But how many eggs should the parasite 
lay and the host accept? For this, we consider a more complex system 
(Fig. 3) where the parasite can lay more than a single egg into a nest 
and where the host can determine a threshold for the number of 
parasitic eggs that will be accepted. Moreover, the host can now 
choose between abandoning the nest or ejecting the parasitic egg, 
which is typical of the considered host species (magpies or 
prothonotary warblers) \ 

We analyse the host-parasite dynamics using large-scale indi- 
vidual based simulations, cf Methods section for more details. We 
observe dynamics reminiscent of cyclic patterns between mafia and 
non-mafia parasites and among rejecter, conditional accepter, and 
accepter hosts. Fig. 4. Mafia parasites are most prevalent when 
acceptance is conditioned upon retaliation and this slowly decreases 
as the host population changes from conditional- accepters to just 
accepters until non-mafia parasites increase in frequency. We also 
observe that almost all mafia parasites return to their hosts. 
Opportunistic or spiteful behaviors emerge only at low frequencies. 
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Figure 4 | Simulation of the complex model of host-parasite interactions depicted in Fig. 3. The dynamics are reminiscent of cyclic pattern between 
mafia and non-mafia parasite and among rejecter, conditional accepter, and accepter hosts. Parameters, G = 1000, Np = 50, = 1000, bh = 5, j^^v = 20, 
Pm = 18, fip = 0.005, fih = 0.01, c„ = 1, = 1, Cp = 1. 



Fig. 4. Opportunistic parasites do not punish, they just return to the 
host, while spiteful parasites punish but do not return. These beha- 
viors take advantage of the increasing accepter hosts in the popu- 
lation, but tend to be unsuccessful since most conditional hosts only 
accept returning parasites who destroy their nest. This suggests that 
all behavioral changes (between the first clutch to the second) require 
retaliation and a second parasitism from the same parasite. 

Furthermore, costs imposed on hosts by the parasite drive changes 
in behavior. In general, nest abandoning occurs when parasites are 
non-retaliators and re-clutching costs are low. Otherwise the host 
just ejects the eggs. Simulations show that as nestling cost increases, 
the frequency of retaliating parasites increases, which is in agreement 
with Eqs. [4]. Parasites have to enforce acceptance with increasing 
nestling cost, since hosts prefer to reject when nestling costs are high. 
Similarly, increasing re- clutching cost increases accepter hosts fre- 
quency, and as a result, mafia parasites decrease in frequency. Fig. 5. 

Interestingly, the majority of the parasites evolved to lay a single 
egg per interaction. Fig. 5, which is in line with observed patterns^ 
This was independent of all costs (parasitism, nestling, and re-clutch- 
ing) with the assumption that there is no extra cost to finding addi- 
tional nests. Intuitively, since a host that accepts / parasitic eggs also 
accepts any lower number, a parasite that lays one egg per nest will 
have the highest acceptance rate. 

Discussion 

Punishing free- riders or individuals who break group rules is 
observed both in animal and human societies^^"^^. However, punish- 
ment by 'rule breakers' or free-riders seems to be less common and 
occurs when individuals use violence or threats to compel another to 
do their bidding. This is most successful when individuals are 
reputed to punish when others refuse to comply^^, for example, mafia 
organizations or school bullies. When the mafia hypothesis was ori- 
ginally proposed by Zahavi in 1979^^ he was opposed on the grounds 
that such behavior is costly for the parasites, without having direct 
benefits^^: mafia behavior would reduce the fitness of rejecting hosts, 
but the parasite's ejected chick would have nothing to gain from the 
parent's action. However, herein, we have shown that retaliatory 
behavior emerges when hosts have plastic behaviors. Thus, while 
the parasite's ejected nestling indeed has nothing to gain from retali- 
ation, the parasite's future nestlings can benefit from this behavior, 
providing an explanation for the retaliation observed in nature^^'^^. 

This explanation may appear surprising: for retaliatory behavior 
to be favored, hosts need to have evolved a plastic response. On the 
other hand, the ability to perform such a plastic behavior is only 
favored if already a large fraction of parasites retaliate. So how can 



both behaviors emerge in the first place? In our model, the host's 
ability to react conditionally comes without exogenous costs and 
benefits. In nature, however, the ability to react conditionally may 
have several additional advantages (e.g. a better defense against nest 
predators), but may also come with additional costs (e.g., as hosts 
need to process surrounding environmental cues). Ultimately the 
emergence of conditional behavior depends on an environment 
where it is beneficial for hosts to build up the capacity to react. As 
our model suggests, once hosts have acquired the ability to react 
conditionally, parasites can use clutch depredation as a way to 
enforce acceptance. 

In experimental studies, retaliatory behavior is hard to distinguish 
from so called farming^ since they both involve depredation as a 
tactic. Whereas mafia parasites use depredation as a response to 
rejection, farming parasites depredate mature clutches before para- 
sitizing them to synchronize the timing with the host, but not to 
induce a specific response. Both create opportunities for the parasite 
by forcing the host to re-clutch. In addition to depredation, retaliators 
monitor the host and ensure acceptance of their egg. As our primary 
interest was in the emergence of mafia behavior and conditional 
responses, we have refrained from including farming behavior (this 
would require a different type of model, that needs to track the age of 
the clutch). However, it is important to note that incorporating farm- 
ing behavior cannot explain host acceptance. Only under retaliation 
can conditional acceptance of the host become beneficial. 

In a previous model, the frequency of accepting hosts depended on 
the success of a second clutch^^. In such a case, the parasite must 
'reward' the accepter hosts by not parasitizing the second clutch. But 
observations show that magpies lay only one clutch during a breed- 
ing season unless the parasite retaliated^^. Our results show that 
accepter hosts will evolve as a response to retaliation and there is 
no need for the parasite to reward accepters, which is in line with 
experimental observations^'^^'^^. As accepter host frequency 
increases, it is no longer necessary for parasites to retaliate. Thus, 
non-retaliators can emerge and increase in frequency. Our model 
also suggests that if retaliation costs are high (eg. reducing total eggs 
laid, small Pm) then the probability that a host is parasitised 
decreases. As a result hosts are more inclined towards rejection 
and thus parasites must retaliate to induce acceptance. Such a situ- 
ation provides a glimpse into the frequency-dependence effects of 
parasitism probability. Hosts will accept if both the frequency of 
retaliatory parasites and the probability of parasitism are high. 
These results corroborate with a previous modeP^, however, we also 
show that the co- evolutionary arms race results in cycles with one 
species lagging behind the other, as shown in Fig. 2. 
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As nestling cost increases the frequency of retaliating parasites increases, (b) As re-clutching cost increases the accepter-hosts increase in frequency, and 
since there is no need to enforce acceptance, mafia parasites decrease in frequency. Parasites again evolved to lay a single egg for all nestling and re- 
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Our results are in line with observed patterns, suggesting that hosts 
learn^'^^'^^ and change their behavior as a response to mafia parasites 
within one breading season^^. It is such learnt behavior from a 
repeated interaction that promotes the evolution of retaliation in 
our model. Typically, repeated interactions and co- evolutionary 
changes are thought to induce parasites to specialize on only one 
or a fewhosts^^. However, our results present retaliation as a behavior 
that the parasite can use to evade specialization and induce accept- 
ance by the host. Other studies which consider manipulative beha- 
vior in parasitic species^^'^^'^^ have suggested a similar effect of 
plasticity^"''*^ A general behavior, such as retaliation, forces a host 
to change it's response and as a result the parasite has a mechanism to 
manipulate multiple hosts without the need to specialise. We further 
speculate that the success of other parasitic species, which can infect 
multiple hosts, may be explained by the same mechanism. 

Methods 

We simulated the evolutionary dynamics of this host-parasite interaction by 
assuming two finite populations of size and Np for hosts and parasites, respectively. 
Our model is based on constant population sizes, which is a good approximation if 
neither hosts nor parasites undergo large fluctuations in their abundance. Relaxing 
this assumption^^'^^'"^^ can lead to a different kind of dynamics which is not considered 
here. In each season a parasite will interact with multiple hosts. All parasites play G 
games repeatedly against randomly selected hosts, resulting in the game payoffs given 
in Fig. 3. In the beginning of each season, the host lays bh eggs in a clutch, which may 
become parasitized. The parasite's behavior is determined by the number of eggs laid 



per host nest bp (with bp < bh), whether it is a retaliator or not, and whether it returns 
to the same host or not. Thus, there are 4:{bh — 1) possible strategies for parasites. The 
host's behavior is determined by the degree of parasitism it can detect or tolerate, t 
(with t<bh), whether mafia parasites are accepted after depredation or not, whether 
returning parasites are accepted or not, and whether the host ejects or abandons the 
entire nest. Fig. 3. For instance, depending on the host's threshold trait t, the host may 
accept the parasitism (if bp < t) incurring a parasitism cost Cp and a nestling cost c„; 
the parasite gains bp eggs. On the other hand, the host can reject the parasitism (if bp > 
t) either by ejecting the eggs or abandoning the nest, either way risking retaliation and 
re-parasitism by the parasite. In the case where these hosts avoid both retaliation and 
re-parasitism, ejector hosts incur Cp while abandoner hosts incur a re-clutching cost, 
Cr. Once rejected the parasite can decide to return, either for retaliation or parasitism, 
and the host can again decide to accept or reject (we assume only ejection in the third 
stage, as unconditional abandoners always have fitness of zero). Overall, this leads to 
Sbh possible strategies for the host. We assume non-overlapping generations in which 
individuals reproduce and contribute offspring to the next generation in proportion 
to their fitness while maintaining constant population sizes for both hosts and 
parasites (this corresponds to a Wright-Fisher process in population genetics). 
Offspring inherit the strategy of the parent. However, at each reproductive step 
mutations can occur with probability fip for the parasite population and fi^ for the 
host population, leading to a random phenotype. In comparison to the minimalistic 
model, where it was not necessary to include mutations, this time the strategies are 
coded by a collection of discrete "genes" that can mutate independently, which 
produces natural transitions from one strategy to the next. 
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